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摘 要 对 于 具有 了 筷 际 结构 和 微 尺度 效应 的 低 渗 透 储 层 ， 其 渗流 现象 的 介 观 机 理 逐 渐 受到 人 们 的 重视 .本 文采 用 微 流 控 
技术 在 硅 片 上 加 工 出 微 通道 作为 微观 模型 , 模拟 岩心 内 部 的 孔隙 结构 , 对 微 通道 中 单 相 和 油水 两 相 的 渗流 规律 进行 了 实验 
研究 . 研究 结果 表明 , 微观 模型 中 带 有 哈 道 段 的 微 通道 存在 启动 压力 梯度 , 这 揭示 了 孔隙 中 的 喉 道 是 低 渗 透 储 层 中 产生 启 
动 压力 梯度 的 主要 介 观 因素 . 启动 压力 与 孔隙 结构 和 驱 替 介质 的 性 质 有 关 , 在 此 基础 上 建立 了 适用 于 单 相 和 油水 两 相 渗 流 
的 启动 压力 预测 模型 。 
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Mesoscopic Mechanism of Threshold Pressure Gradient in Low 

Permeability pe 
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(State Key Laboratory of Multiphase Flow in mt ing, Xi’an Jiaotong University, Xi’an 710049, China) 
S 


Abstract The mesoscopic mechanism of flow phenomena in low permeability reservoirs has 


been paid considerable attention due complex pore structure and the micro-scale effect. In 


this paper, micro-channels are pr in silicon chips by microfluidic technology to simulate the 
pore structure inside low perpoRot reservoirs. The experiments are conducted on the porous flow 
performance of single-phas bil-water two-phase flows by using the micro-channels manufactured. 
The results show that Re pressure gradient appears in the micro-channels with throats, 
which indicates that ee in a micro-channel is the main nn resulting in threshold 
pressure gradient. The threshold pressure is associated with the structire of pore and the property of 
displacing medium. On this basis, a new model is proposed frredicting the threshold pressure of 
single-phase and oil-water two-phase flows in low permeabi porous media. 


Key words threshold pressure gradient; ro LO ter two-phase flow; microchannel 
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XX 
在 我 国 探 明 的 原油 储量 中 ， 低 渗透 油田 电 吝 间 
常 大 的 比例 , 是 陆地 石油 工业 持续 发 展 的 重要 储备 。 
启动 压力 梯度 对 油 藏 采 出 程度 具有 显著 影响 四， 因 
此 对 低 渗透 油 藏 启动 压力 梯度 机 理 的 研究 意义 重大 . 
由 于 低 渗 透 储 层 一 般 具有 了 筷 喉 细小、 比 表面 积 大 、 贾 
敏 效应 和 表面 分 子 力作 用 强烈 等 特点 ， 导 致 难以 精 
确 地 描述 其 流动 机 理 目 . 为 此 ,国内 外 学 者 进行 了 
大 量 的 研究 . Sun 等 里 对 岩石 、 玻 璃 微观 模型 的 多 
相 渗流 研究 表明 ， 低 渗透 油 藏 的 渗流 规律 、 远 不 同 
于 高 渗透 油 藏 . 宋 付 权 等 由 对 油水 两 相 流 的 研究 表 
明 , 在 低 渗 透 油 藏 中 , 启动 压力 梯度 不 会 对 油 相 相 对 
渗透 率 产生 影响 ,但 会 对 水 相 相 对 渗透 率 产生 不 同 
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程度 的 影响 . 黄 延 章 等 回 的 实验 研究 表明 ， 低 渗透 
油田 中 启动 压力 梯度 普遍 存在 , 且 岩 心 渗透 率 越 低 ， 
其 启动 压力 梯度 越 大 ; 含水 饱和 度 越 小 , 启动 压力 梯 
度 越 大 . Thomas [9 探讨 了 启动 压力 与 岩心 渗透 率 、 
迁 曲 度 、 和 孔隙 度 、 形 态 因 子 和 流体 表面 张力 的 关系 ， 
并 回归 得 到 了 经 验 公式 , Ito 等 中 分 别 采用 逐步 递增 
法 、 余 压 法 、 常 规 压 汞 法 ,测定 了 兰 样 的 启动 压力 ， 
结果 表明 , 余 压 法 获得 的 启动 压力 值 小 于 逐步 递增 
法 的 测定 值 . Huet 加 基于 寡 律 模型 ， 采 用 非 线性 最 
小 二 乘法 拟 合 得 到 岩心 渗透 率 、 启 动 压力 与 岩 样 物 
性 特征 的 相互 关系 , 并 结合 压 汞 法 得 到 的 实验 数据 ， 
分 析 得 到 了 较为 符合 实际 的 经 验 公 式 . 刘 丽 等 四 通 
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过 实验 和 理论 分 析 相 结合 的 方式 建立 了 低 渗 透 非 均 
质 砂 岩 油 藏 启动 压力 梯度 确定 方法 ， 能 够 用 于 有 效 
指导 低 渗 透 非 均 质 砂 宕 油 藏 的 合理 井 距 的 确定 。 

岩心 的 渗流 实验 是 低 渗透 油 藏 开发 中 最 早 使 用 
也 是 目前 最 常 使 用 的 技术 ， 它 反映 了 流体 在 多 孔 介 
质 中 流动 的 宏观 统计 平均 效果 ,无 法 真实 地 反映 其 
内 部 的 细 观 渗流 特性 . 微 蚀 刻 技 术 的 发 展 为 细 观 渗 
流 特 性 的 研究 提供 了 技术 平台 ， 利 用 该 技术 刻 蚀 出 
含有 孔道 - 喉 道 的 硅 - 玻璃 微观 模型 来 模拟 真实 
的 岩心 样本 ， 并 与 数字 显 微 摄像 技术 和 微 流体 测试 
技术 相 结 合 ， 可 以 达到 对 渗流 过 程 的 定量 研究 ， 揭 
示 出 孔 除 - 喉 道 结构 对 低 渗透 介质 启动 压力 梯度 的 
影响 ， 


1 模型 建立 
根据 Purcell L109] 模型 Huet [9 通过 分 析 研 究 启 


动 压力 的 双 曲 线 特征 ， 建 立 了 与 岩心 物性 参数 、 启 
动 压力 、 润 湿 相 饱和 度 相 关 的 渗透 率 估算 模型 即 : 


1.65 
= 1017003 二 大 (ze) 0 SR 


P. = 640.05%0.82 开 一 82( 
kPa; 入 为 


入 一 0.00989%-0.63 开 037( 
式 中 , K 为 渗透 率 , mDi 


结构 系数 ;Sw wate 爱 ; 4 为 岩心 孔隙 度 . 
Thomas 1 从 毛细 管 模型 、 泊 肃 叶 流 模型 出 发 
建立 了 启动 压力 的 估算 模型 ， 
IC 1 /9 
及 = 10075- 记 7V 天 (4) 
在 此 基础 上 , 我 们 建立 了 如 下 模型 : o: 
SS 
P, 0.44r-0.5736 二 b 3.9709 yo 0 0.5 GN 
(5) 


式 (4)、(5) 中 , PP 为 启动 压力 , kPa; o 为 流体 表面 张 

力 , mN/m; 7 为 孔 际 形状 因子 ;7 为 岩心 迁 曲 度 ; 9 为 
岩心 孔 际 度 ; K 为 渗透 率 , Di jw 为 黏度 , mPas; Gw 
为 含水 率 ， 


2 启动 压力 测试 


2.1 实验 装置 

0 1 所 示 的 高 压 实验 装置 上 完成 . 实 
验 台 以 高 纯 压 缩 氮气 作 动力 ， 推 动 储 液 饶 内 流体 进 
pd 采用 压 差 - 流量 法 求 
得 其 启动 压力 . 实验 以 模拟 油 、 日 油 、 藻 馏 水 和 油水 
两 相 流 (6B 为 白 油 - 水 两 相 流 中 的 含水 率 ) 作 驱 蔡 
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介质 ， 其 物性 参数 如 表 1 所 示 。 所 用 的 岩心 为 延长 
油田 不 同 区 块 的 天 然 岩 心 , 其 物性 参数 见 表 2 所 示 .， 
如 图 2(a) 所 示 为 实验 用 微观 模型 中 的 孔道 - 喉 道 结 


构 , 根据 真实 岩心 的 空 孙 结 构 由 化 学 刻 蚀 技术 得 到 ， 
其 几何 尺寸 测量 结果 如 表 3 所 示 ， 其 中 前 7 条 通道 
都 带 有 一 个 喉 道 段 , 最 后 
没有 设置 喉 道 段 。 微观 模型 中 的 孔道 - 
互 并 联 , 如 图 2(b) 所 示 。 


一 条 通道 T-4-1 为 直通 道 ， 
喉 道 结构 相 


| 热电 偶 
E 风 | 
| 稚 上 网 。 压力 传感器 
所 过 滤器 
气 HH 全 = 
并 守信 咯 (证 | -NA 


流 | 岩心 夹 持 器 阳 了 大 
图 1 实验 台 系 统 , (a) 微 通 道 渗流 实验 台 (b) 岩心 渗流 实验 台 


Fig. 1 Experiment@ / S 
microchanngls (b) Test devices for porous flow in cores 


(a) Test devices for porous flow in 


\ 
x 表 1 流体 工 质 的 物性 参数 


raf 1 Physical properties of fluid working fluid 
A 物质 / 密度 / 表面 张力 / 


所 

kg:m-3 mPa:s mN-.m-! 
998.3 0.987 71.80 
ss 833.2 9.76 29.10 
模拟 油 785.1 2.33 25.40 
Bw =90% 980.2 1.25 47.32 
Pw =80% 961.2 1.18 52.57 
Bw=70% 943.2 1.30 63.10 
Bw=60% 920.5 1.45 56.45 
Bw=50% 901.3 1.66 53.39 
Bw=40% 880.0 49.8 36.72 
Bw=30% 860.2 36.2 33.16 
Bw=20% 840.5 12.5 30.05 
Bw=10% 821.2 11.0 28.99 


2.2 实验 步骤 

微观 模型 渗流 实验 步骤 如 下 : 1) 采用 恒温 烘 干 
法 干燥 模型 ，2) 将 模型 抽 真 空 至 驱 替 工 质 饱和 ; 
3) 微观 模型 连 入 试验 段 管 路 ， 接 通 前 ， 排 出 装置 中 
残留 的 空气 ; 4) 采用 稳定 的 驱 蔡 压力 ， 驱 蔡 数 小 时 
至 其 形成 稳 态 渗流 ， 用 电子 天 平 收 集 驱 蔡 出 的 流体 
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并 记录 时 间 ; 5) 调整 压力 ,并 重复 步骤 小 。 岩心 渗 
流 过 程 类 似 ， 


表 2 低 渗 透 岩心 物性 参数 
Table 2 Physical properties of low permeability 


cores 

岩心 长 度 / 直径 / 孔隙 度 / 渗透 率 / 
编号 cm cm % 10-3um2 
1## 4.432 2.532 12 1.17 
2# 4.422 2.532 13.1 3.26 
3# 4.412 2.532 14.2 5.1 
4# 4.958 2.532 14.8 5.23 
5# 5.672 2.532 15.2 11.5 
6# 5.524 2.532 15.5 12.7 
7?# 5.312 2.532 20.2 38.9 
8## 4.774 2.532 20.9 56.6 
9# 5.356 2.532 12.5 1.23 
10## 4.956 2.532 12.9 3.31 
11## 5.326 2.532 

12# 4.934 2.532 

13## 5.629 2.532 

14# 6.392 2:532 

15# 4.996 2.532 

16# 5.329 2.532 


通道 并 联结 构 
tructure of parallel 
channels 


(a) 孔道 - 喉 道 单 通道 结构 

(a) Structure of a channel NS 
with a throat 

图 2 微观 结构 模型 


Fig. 2 Microstructure model 


表 3 微 通道 几何 尺寸 


Table 3 Geometries of microchannels 


通道 。 孔道 宽度 / 喉 道 宽度 / 孔 喉 比 叭 道 长 亡 / 
编号 hm hm 

T-1-1 195 31 6.3 

T-1-2 195 41 4.8 96 
T-1-3 195 51 3.8 98 
T-2-1 192 37 5.1 101 
T-2-2 192 45 4.2 99 
T-2-3 192 57 3.4 102 


T-3-1 167 21 7.9 104 
T-4-1 190 一 


3 实验 结果 与 分 析 
图 3 为 T-4-1 微 通道 流体 速度 与 压力 梯度 的 关 
系 , 其 中 微 通道 T-4-1 为 直通 道 , 其 上 没有 设置 喉 道 
段 . 从 图 3 中 可 以 看 出 , 流体 流速 与 压力 梯度 呈正 
比例 关系 , 符合 达 西 定律 , 不 存在 启动 压力 梯度 . 
图 4 为 微观 模型 中 单 相 水 渗流 的 启动 压力 梯度 


与 渗透 率 关 系 曲线 ,其 中 微观 模型 是 指 表 3 中 给 出 
的 前 7 条 微 通道 所 在 的 硅 - 玻璃 微观 模型 ,其 上 痢 
分 别 分 布 有 一 个 喉 道 段 . 从 中 可 以 看 出 , 启动 压力 梯 
度 与 渗透 率 呈 线性 关系 , 显然 不 同 于 岩心 模型 分 析 
原因 ， 岩心 模型 中 孔隙 - 喉 道 结构 的 数量 很 多 而 且 
分 布 比 较 复杂 ， 而 微观 模型 中 只 有 一 个 喉 道 段 ， 而 
且 不 存在 孔隙 - 喉 道 结构 分 布 不 确定 的 问题 ， 其 流 
动 阻力 主要 来 自主 要 是 孔 除 和 喉 道 之 间 变 化 引起 的 
局 部 阻力 以 及 在 狭 窗 的 喉 道 段 中 的 流动 阻力 . 但 是 ， 
在 两 种 模型 中 启动 压力 梯度 随 渗透 率 增加 而 减 小 的 
趋势 是 一 致 的 . 与 图 3 所 作出 的 分 析 对 比 得 知 ， 岩 
心 模型 中 存在 启动 压力 梯度 是 由 于 其 中 存在 和 孔 际 - 
喉 道 结构 而 造成 的 。 


3 


AP.Z-VMPa.mr-l 
图 3 T- 人 4 水 通道 流速 与 压力 梯度 的 关系 


Fig. 3 Relationship between flow velocity and pressure 
C gradient in channel T-4-1 


\ 


0 100 200 300 400 500 600 
K/D 


图 4 微观 结构 模型 的 启动 压力 梯度 与 渗透 率 关 系 曲线 


Fig. 4 Relationship between threshold pressure gradient and 


permeability of microscopic models 


根据 Thomas 外 提出 的 渗透 率 估算 模型 ,通过 
分 析 研 究 , 我 们 发 现 启动 压力 还 与 驱 替 流体 饱和 度 、 
黏度 等 参数 有 关系 。 在 此 我 们 采用 寡 律 模型 对 其 进 
行 拟 合 ,得 到 新 的 启动 压力 预测 模型 ， 见 式 (5)。 如 
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图 5 所 示 ， 当 启动 压力 较 大 时 ， 实 验 值 与 拟 合 值 吻 
合 得 很 好 , 拟 合 曲线 的 相关 系数 R? 在 0.98 以 上 ; 启 
动 压力 较 小 时 ， 计 算 值 远 小 于 实验 值 ， 这 很 可 能 是 
由 于 油水 两 相 流 含水 率 较 大 、 黏 度 较 小 时 二 者 的 影 
响 减 小 而 预测 模型 没有 考虑 这 一 点 造成 的 。 因 此 可 
以 认为 拟 合 结果 是 可 靠 的 。 

图 6 给 出 了 传 水 率 为 10%~90% 时 微观 模型 启动 
压力 的 试验 值 与 计算 值 (以 通道 T-3-1 为 例 ). 从 比较 
结果 来 看 , 在 含水 率 为 10%~60% 的 条 件 下 , 渗流 的 
启动 压力 较 高 , 计算 值 与 实验 值 有 较 好 的 相关 性 , 而 
含水 率 为 70%~90% 时 ,渗流 的 启动 压力 低 , 计算 值 
远 小 于 试验 值 ， 这 很 可 能 是 由 于 启动 压力 不 存在 造 
成 的 。 这 表明 本 文 提出 的 启动 压力 预测 模型 ， 在 一 
定 范围 内 能 够 有 效 预 测 低 渗透 储 层 的 启动 压力 . 


油水 两 相 流 


Pc/kPa 


PN 
图 5 微观 模型 启动 压力 的 预测 


Fig. 5 Prediction on threshold pressure of a microscopic 


model 


全 油水 两 相 流 


一 人 一 理论 曲线 


Poe/kPa 


0 人 人 
10 20 30 40 50 60 70 80 90 
GB,/% 


图 6 微观 模型 的 油水 两 相 流 启动 压力 与 含水 率 的 关系 


Fig. 6 Relationship between threshold pressure and moisture 


content of oil-water two-phase flow in a microscopic model 


图 7 为 岩心 样品 实测 的 启动 压力 与 理论 模型 计 
算 值 的 对 比 图 . 结合 岩心 样品 实测 的 启动 压力 和 渗 
透 率 , 并 与 Thomas 向 、Huet [6 理论 预测 式 进 行 对 


比分 析 ， 如 图 7 (a)~(c) 所 示 。 结果 表 明 ，Huet 模 
型 可 以 精确 地 预测 储 层 渗透 率 ， 其 相对 误差 不 大 于 
土 10%。 从 图 7 (a) 可 以 看 出 ， 当 岩心 渗透 率 较 低 时 ， 
水 、 模 拟 油 在 岩心 中 的 启动 压力 相差 较 大 ; 当 兰 心 的 
渗透 率 较 高 时 , 它们 的 启动 压力 基本 一 致 , 可 知 表面 
张力 等 微观 特性 在 低 渗 透 率 储 层 的 影响 程度 远大 于 
高 渗透 率 储 层 .Huet 和 Thomas 模型 的 启动 压力 预 
测 结果 与 实验 值 误差 较 大 ， 因 此 ， 上 述 两 种 模型 都 


实验 值 


田 水 


OO 模拟 油 


线 
s 预 测 曲线 
omas 预 测 曲 线 


(a) 启动 压力 预测 模型 
(a) Threshold pressure predictionmodel 


四 水 
® pn 


10 对 误差 区 间 Ee 
有 
> KK 


1 YL65 
由 442 


2 Kg=1017003 


(IS) Gt 


30 
K./mD 
(b) 岩心 渗透 率 预 测 模 型 


(b) Permeability prediction modelof cores 


50 60 


四 水 
© 模拟 油 


Pe-w=Pc-ex 
cp 1 
Pe m=0.440 0%5736 ato he 


5 
时 _05,3 -1.5962 
K-03614 -0513 DD, 


» 
100 150 200 250 300 350 400 
Pca /kPa 
(c) 启动 压力 实测 值 与 计算 值 的 对 比 
(c) Comparison of measured andcalculated values of 
threshold pressure 


0 350 


ChinaXiv 合 


4 期 刘 冬 旭 等 : 低 渗透 介质 启动 压力 梯度 的 介 观 机 理 研究 


4) 对 于 相同 渗透 率 的 岩心 模型 , 渔 湿 相 的 启动 


es 压力 梯度 略 高 于 非 润 湿 相 ， 并 随 渗 透 率 的 降低 而 升 


高 ， 两 者 成 反比 例 函 数 关系 ; 在 带 有 了 筷 际 - 喉 道 结 
构 的 微观 模型 中 ,启动 压 力 梯 度 随 渗透 率 的 降低 而 
升 高 ， 且 两 者 成 线性 关系 . 带 有 了 筷 际 - 喉 道 结构 的 
微观 模型 能 够 在 一 定 程度 上 说 明 岩 心 模型 中 孔 际 - 
喉 道 结构 对 启动 压力 梯度 的 影响 ， 
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4 结 论 


1) 孔隙 - 喉 道 微观 模型 中 存在 着 启动 压力 ， 基 
于 Thomas 模型 得 到 的 新 模型 拟 合 关系 预测 出 的 启 
动 压力 更 接近 于 实测 数据 ,准确 度 较 高 。 

2) 微观 模型 中 ,与 没有 喉 道 段 的 直通 道 相 比 带 
有 喉 道 段 的 通道 存在 启动 压力 梯度 ， 而 且 其 中 喉 道 
尺寸 的 减 小 ， 会 导致 启动 压力 梯度 急剧 增加 ， 这 很 
好 地 证 明了 和 扎 阶 中 狭 罕 的 喉 道 是 启动 压力 梯度 存在 
的 主要 因素 . 

3) 低 渗 透 油 藏 中 存在 着 启动 压力 。 本文 提出 的 
启动 压力 预测 模型 ， 实 现 了 对 启动 压力 较为 准确 地 
预测 ， 当 岩心 绝对 渗透 率 低 于 一 定数 值 时 ，Huet 模 
型 可 以 精确 地 预测 储 层 的 渗透 率 ， 
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